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Abstract 

We discuss the bounds on the cosmological lepton number from Big Bang 
Nucleosynthesis (BBN), in light of recent evidences for a large value of the 
neutrino mixing angle ^13, sin^ ^^13 > 0.01 at 2-a. The largest asymmetries for 
electron and /i, r neutrinos compatible with ^He and primordial yields are 
computed versus the neutrino mass hierarchy and mixing angles. The flavour 
oscillation dynamics is traced till the beginning of BBN and neutrino distri- 
butions after decoupling are numerically computed. The latter contains in 
general, non thermal distortion due to the onset of flavour oscillations driven 
by solar squared mass difference in the temperature range where neutrino 
scatterings become inefficient to enforce thermodynamical equilibrium. De- 
pending on the value of ^13, this translates into a larger value for the effective 
number of neutrinos, A^eff- Upper bounds on this parameter are discussed for 
both ncTitrino mass hierarchies. Values for Ncs which are large enough to 
be detectable by the Planck experiment are found only for the (presently 
disfavoured) range sin^ ^13 < 0.01. 
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1. Introduction 

Nowadays flavour neutrino oscillations are well established thanks to a 
plethora of experimental results on the detection of reactor, accelerator, at- 
mospheric and solar neutrinos. Two neutrino mass-squared differences and 
three mixing angles drive oscillations among the three active neutrinos, all of 
them, except one angle, known with a precision better than 25% PQ121E]- The 
last parameter which is not so well known is the mixing angle ^13. Until very 
recently only an upper bound on sin^ 6'i3 existed, while the present year we 
have witnessed the first indications of non-zero 6*13 values from the analysis of 
global data |2], |3] , especially recent z/^ — )■ searches at the T2K long-baseline 
experiment This result opens up the possibility of measuring, in a not 
so distant future, the pattern of neutrino masses (the mass hierarchy) and a 
possible CP violation in the leptonic sector [5j , the last remaining unknowns 
together with the nature of neutrinos (Dirac or Majorana). 

Neutrino oscillations have implications in many research areas in particle 
and astroparticle physics. However, the consequences of non-zero neutrino 
mixing in cosmology are not so important, despite the fact that relic neu- 
trinos are the second most abundant particles in the Universe, with almost 
the same number density as photons. The reason is well known: in first 
approximation all neutrino flavours were produced by frequent interactions 
in the early hot Universe, with the same momentum spectra. Thus neutrino 
oscillations, although effective right after neutrino decoupling, do not modify 
the properties of cosmological neutrinos (except for very small effects from 
non-instantaneous neutrino decoupling p]). This holds for active neutrino 
oscillations, while the case of active-sterile oscillations does lead to effects on 
cosmological observables if one or more extra neutrino species are populated 
(see e.g. [3 El El [ID]). 

There exists, however, one situation where active-active oscillations have 
an impact on cosmological neutrinos, namely when a large flavour neutrino 
asymmetry was previously created. It is usually assumed that such an asym- 
metry, parameterized by the number density ratios 

= , a = e, /i, r , (1) 

should be of the same order of the cosmological baryon number r/^ = (r?,;, — 
ni)/n^, due to the equilibration by sphalerons of lepton and baryon asymme- 
tries in the very early universe. Thus one does not expect flavour neutrino 
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asymmetries much larger than a few times 10~^°, the value of rjh measured by 
present observations, such as 7-year data from the WMAP satellite and other 
cosmological measurements [TT]. There are, however, some models where a 
lepton asymmetry orders of magnitude larger than the baryon one could sur- 
vive, see e.g. [121 US], in the neutrino sector with an influence on fundamental 
physics in the early universe, such as the QCD transition [H] or a potential 
relation with the cosmological magnetic fields at large scales [T5] . 

Present cosmological observations are not sensitive to a large neutrino 
asymmetry if {rji^l < 10~^. Only larger values lead to a significant enhance- 
ment of the contribution of active neutrinos to the radiation energy density 
of the Universe or to changes in the production of light elements in Big Bang 
Nucleosynthesis (BBN). In particular, the primordial abundance of ^He de- 
pends on the presence of an electron neutrino asymmetry and sets a stringent 
bound on rji,^ which does not apply to the other flavours unless neutrino os- 
cillations are effective before BBN, leaving a total neutrino asymmetry of 
order unity unconstrained [Ml [IT] . 

A decade ago, it was shown that flavour neutrino conversions in the early 
Universe are indeed suppressed by matter effects at large temperatures, and 
it is only at temperatures T < 10 MeV that oscillations set on and are 
large enough to achieve strong flavor conversions before BBN [TH [191 I20l 
1^ . For the present measured values of neutrino mixing parameters, the 
degree of flavour equilibration depends on the value of 613. As discussed 
in [22], this parameter fixes the onset of flavour oscillations involving z/g's, 
which in turn determines whether neutrinos interact enough with electrons 
and positrons to transfer the excess of energy density due to the initial r]^^ 
to the electromagnetic plasma. Recently we have found the BBN bounds 
on the cosmological lepton number for a range of initial flavour neutrino 
asymmetries [23]. 

Prompted by the recent indication of non-zero values for ^13 and the hints 
of possible extra radiation from cosmological data [11] , we have updated the 
analysis in [23] with the aim of finding the BBN bounds on both the total 
neutrino asymmetry and its maximum contribution to the radiation content 
of the Universe in the whole range of 613 values allowed by oscillation data, 
as well as considering both neutrino mass hierarchies. 

This letter is organized as follows. We introduce in Section |2] the formal- 
ism of kinetic equations which rule the evolution of neutrino distributions 
and describe the dynamics of neutrino asymmetries in the epoch just imme- 
diately BBN. In Section [3] we study the BBN constraints on lepton number 
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that can obtained for a wide range of initial neutrino asymmetries, with 
emphasis on the effects of a mixing angle ^13 in the region experimentally 
favoured. Finally in Section |4] we give our concluding remarks. 

2. Neutrino evolution in the presence of lepton asymmetries 

Active neutrinos were produced in the very early Universe and their en- 
ergy spectrum is kept in chemical and kinetic equilibrium by weak inter- 
actions until temperatures T ~ O(MeV), when the corresponding collision 
rates fall below the cosmological expansion rate. Therefore, if flavour neu- 
trino asymmetries existed, at larger temperatures the neutrino distribution 
of momenta is a Fermi-Dirac spectrum parameterized by a temperature T 
(the same of e~^e~ and photons) and a well defined chemical potential fi,^^ 
for a = e, yU, r. Each flavor neutrino asymmetry in Eq. Q can be expressed 
in terms of the corresponding degeneracy parameter = l^v^/T as 

= ^2^(^'ea + e') , (2) 

with C(3) ^ 1.20206. This expression is modified later by a factor (T^^/T^)^, 
when e+e^ pairs annihilate into photons. The corresponding contribution of 
neutrinos in equilibrium to the total energy density, usually parameterized as 
Pr/p-y = 1 + 7 /8{A/liy^^Ncs after the e~^e~ annihilation phase, is enhanced 
for non-zero neutrino asymmetries as follows 

iVeff = 3+ J2 

The parameter Nes is the effective number of neutrinos whose standard value 
is 3 in the limit of instantaneous neutrino decoupling. 

We are interested in calculating the evolution of the active neutrino spec- 
tra from large temperatures, when they followed a Fermi-Dirac form, until 
the BBN epoch. This includes taking into account neutrino interactions 
among themselves and with charged leptons, as well as flavor oscillations, 
which become effective at similar temperatures. In such a case the best way 
to describe neutrino distributions is to use matrices in flavor space Qp [211 125] . 
For three active neutrino species, we need 3x3 matrices in flavor space gp 
for each neutrino momentum p, where the diagonal elements are the usual 
occupation numbers and the off-diagonal ones encode phase information and 
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vanish for zero mixing. The corresponding equations of motion (EOMs) for 
are the same as those considered in references [221 ES] , where the reader 
can find more details on the approximations made to solve them and related 
references, 



and similar for the antineutrino matrices Pp. The last term corresponds to 
the effect of neutrino collisions, i.e. interactions with exchange of momenta, 
which are implemented as in ref. [22] • These collision terms, proportional to 
the square of the Fermi constant Gp, are crucial for modifying the neutrino 
distributions to achieve equilibrium with and, indirectly, with photons. In 
the absence of neutrino mixing, the EOMs include only collision terms and 
preserve the flavour neutrino asymmetries r]^^. 

The first term on the right-hand side of Eq. ^ describes flavor oscilla- 
tions. 



where p = |p| and M is the neutrino mass matrix (opposite sign for antineu- 
trinos), which in the flavour basis is not diagonal and includes the mixing 
parameters that characterizes the vacuum term of oscillations. In our calcu- 
lations we have fixed both mass-squared differences and the angles 612 and 
623 to the best-fit values in [3] . Varying these parameters within the allowed 
3(T ranges does not modify our results. Instead, we will consider the whole 
presently allowed range of 613 values, approximately from 0.001 to 0.035—0.05 
at 3(7 (depending on the reactor neutrino fluxes, see [21 [3]), adding the case 
of zero 613 for comparison. 

Matter effects are included via the term proportional to the Fermi con- 
stant Gp, the so-called neutrino potentials. The one proportional to g — 
g, where g = J Ppd^p/(27r)^ (and similar for antineutrinos) arises from 
neutrino-neutrino interactions and it was shown in [THl UHl EQ] that for the 
relevant values of neutrino asymmetries this matter term dominates but does 
not suppress flavour oscillations. Instead, this term leads to synchronized os- 
cillations of neutrinos and antineutrinos of different momenta [2U] . 

The small baryon density in the early Universe implies that the refractive 
matter term proportional to the difference between the charged lepton and 
antilepton number densities (CP asymmetric) can be neglected compared to 
the CP symmetric term, proportional to the sum of energy densities [2ll |27] . 



. dgp 



Pp, ^'p] + C[gp, gp] 



(4) 



dt 
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This appears in Eq. Q with E, the 3x3 flavor diagonal matrix of charged- 
lepton energy densities. The onset of flavour oscillations occurs at a the 
temperature at which the vacuum and charged-leptons background terms 
become equal in magnitude, i.e. when the following terms 

Am' 8V2Gfp 

are comparable for the relevant mixing parameters Am^ and 6, and the 
charged- lepton energy density pi (see e.g. [IS]). Tau leptons are too heavy 
to have a significant density at MeV temperatures, while the energy density 
of fi^ is exponentially suppressed, leading to — mixing driven by Amg^ 
and ^23 at T ~ 15 MeV IlISJ, when weak interactions are fully effective. Thus, 
as in previous works, our numerical calculations start at T = 10 MeV with 
initial degeneracy parameters = (.^ = and ^e- 

For flavour neutrino oscillations involving z/g's, the crucial parameters are 
Amgj^ and ^13 because they fix the moment when the neutrino conversions 
become important, i.e. when the absolute values of both terms in Eq. ^ are 
equal. One easily finds that this occurs at a temperature 
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T, ^ 19.9 (^^j (^^^ j MeV (7) 

for cos 2^13 ~ 1 and taken as relativistic particles. For |Am|^| = 2.5 x 10~^ 
eV^ and an average neutrino momentum, one finds ~ 5 MeV. If Amg^^ > 
(normal neutrino mass hierarchy, NH) both terms in Eq. (|6| have the same 
sign and neutrino oscillations follow an MSW conversion when the vacuum 
term overcomes the matter potential at T ~ Tc. The degree of conversion 
depends in this case on the value of 613 [HI [191 EO], being very efficient 
compared with ^13 = if this mixing angle presents a value close to the 
upper bound, as can be seen in figure [T] for one particular case with zero 
total lepton number and different choices of 613. 

The conversion for non-zero ^13 is more evident for the inverted mass 
hierarchy, as shown in figure [T| due to the resonant character of the MSW 
transition for Am'^^ < 0. Indeed, for IH the sum of the two terms in Eq. 
^ vanishes and the equipartition of the total lepton asymmetry among the 
three neutrino flavours is quickly achieved, even for sin^^i3 < 0.01, unless of 
course 6*13 is extremely small. Finally, for negligible ^^13 flavour oscillations 



6 



0.6 I , , , r 




-1 I ' ' ' ' ' ' ' ' 1 

10 7 5 3 2 1 

(MeV) 

Figure 1: Evolution of the flavor neutrino asymmetries when r/Jf = —0.82 and zero total 
asymmetry. The outer solid curves correspond to vanishing ^13 (black lines), while the 
inner ones (red lines) were calculated in the NH for two values of sin^^ia: from left to 
right, 0.04 and 0.02. The same two values of sin^ ^13 apply to the cases shown as blue 
dotted lines, but in the IH. 



are not effective until T < 3 MeV (outer lines in figure [T|, a value that can 
be found substituting Aml^ for Aml^ = 7.6 x 10~^ eV^ in Eq. 

The moment when flavour oscillations in the presence of neutrino asym- 
metries become effective is important not only to establish the electron neu- 
trino asymmetry at the onset of BBN, but also to determine whether weak 
interactions with e^e~ can still keep neutrinos in good thermal contact with 
the ambient plasma. Oscillations redistribute the asymmetries among the 
flavours, but only if they occur early enough would interactions preserve 
Fermi-Dirac spectra for neutrinos, in such a way that a chemical potential 
/ij,^ is well defined for each i]^^ and the relations in Eqs. ^ and ^ remain 
valid. For instance, if the initial values of the flavour asymmetries have oppo- 
site signs, neutrino conversions will tend to reduce the asymmetries which in 
turn will decrease iVefj. But if flavour oscillations take place at temperatures 
close to neutrino decoupling, this would not hold and an extra contribution 
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Figure 2: Final contribution of neutrinos to the total radiation energy density, 
parametrized with N^s, as a function of the total neutrino asymmetry for a particular 
value of the initial electron neutrino asymmetry (?7jf — —0.82). From top to bottom, the 
various lines correspond, respectively, to the following cases: no neutrino oscillations (?7i/^ 
conserved), 6*13 — 0, and sin^ — 0.04 for normal (red solid line) and inverted (blue 
dotted line) neutrino mass hierarchy. 

of neutrinos to radiation is expected with respect to the value in Eq. ([s]), as 
emphasized in |22] . 

A way to see the role of flavour oscillations on the reduction of the final 
value of Ncs from neutrino asymmetries is given in figure [2] Here we have 
fixed the initial electron neutrino asymmetry to 77™ = —0.82 as in figure [l| 
but varied the total asymmetry in the range —0.8 < r]u ^ 0.8. In the absence 
of neutrino mixing the final value of N^s is that given by Eq. ([s]), directly 
related to the chemical potentials, and for this particular range it can be 
as large as A^efi — 4.6. Instead, when oscillations are included the three 
flavour asymmetries are modified and the contribution of neutrinos is largely 
reduced, even for = 0. Finally, for sin^ ^13 = 0.04 and both NH or IH, the 
final flavour asymmetries are given by rj^^ ~ rji^/S. In such a case, we expect 
neutrinos to almost follow Fermi-Dirac spectra and N^ff as given in Eq. (|3|. 
For instance, for r/,^ = 0.8 one expects C,ue^,T — 0-38 and a total contribution 
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to the radiation energy density of N^g ^ 3.19, very close to what we find in 
our numerical calculations for sin^6'i3 = 0.04, while we found Ne^ ~ 3.43 for 

^13 = 0. 

3. BBN bounds on the cosmological lepton asymmetry 

Cosmological neutrinos influence the primordial production of light el- 
ements in two ways. First of all, they contribute to the radiation energy 
density that fixes the expansion of the Universe during BBN, a background 
effect that is parameterized with N^s- For the particular case of ^He, its 
primordial abundance is enhanced for N^s > 3, as in the case of neutrino 
asymmetries. On the other hand, electron neutrinos and antineutrinos are 
involved in the charged current weak processes which rule the neutron/proton 
chemical equilibrium. Thus any change in the u^, or u^. momentum distribu- 
tion can shift the neutron/proton ratio freeze out temperature and in turn 
modify the primordial ^He abundance. This is what happens for a non-zero 
I'e — i^e asymmetry at BBN, that shifts the neutron fraction towards larger 
or smaller values for negative or positive values of rj^^, respectively. 

We have performed an analysis of the effects of flavour neutrino asym- 
metries on the BBN outcome in a similar way as in [23j. We have solved 
the EOMs described in the previous section in a wide range of values for the 
total lepton (neutrino) asymmetry r^^, unchanged by oscillations, and the 
initial electron neutrino asymmetry rj™. The obtained time-dependent neu- 
trino distributions are then used as an input for the pubhc numerical code 
PArthENoPE [2HI 122] , as described in |23] • The computed abundances of both 
the ratio ^H/H and the ^He mass fraction, Yp, are compared with the corre- 
sponding experimental determinations in order to find the allowed region of 
asymmetries ri,y and Here, as in [23], we consider for the primordial 
abundance the value 

^H/H = (2.87 ± 0.22) x 10"^ , (8) 

obtained by averaging seven determinations from different Quasar Absorp- 
tion Systems ^U\. For the ^He mass fraction we use, as a conservative value, 
the result of the data collection analysis performed in [30] , 

Yp = 0.250 ± 0.003 . (9) 

The main results of our BBN analysis are shown in figure [3]for the adopted 
determinations of and ^He. From this plot one can easily see the effect of 
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Figure 3: 95% C.L. contours from our BBN analysis in the iji, — 77™ plane for several 
values of sin^6'i3: (black solid line), 0.04 and NH (red solid line), 0.04 and IH (blue 
dotted line). The case of no neutrino flavour oscillations is shown for comparison as the 
black dashed contour. 
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Figure 4: The shadowed region corresponds to the values of the total neutrino asymmetry 
compatible with BBN at 95% C.L., as a function of dy^ and the neutrino mass hierarchy. 

flavour oscillations on the BBN constraints on the total neutrino asymmetry. 
In the absence of neutrino mixing the value of rj,^^ is severely constrained by 
^He data, arising from a narrow region for the electron neutrino degeneracy, 
—0.018 < < 0.008 at 68% CL. Instead, the asymmetry for other neutrino 
flavours could be much larger, since the absolute value of total asymmetry 
is only restricted to the region \ri^\ < 2.6 [161 [IZ]- As we have previously 
seen, flavour oscillations modify this picture and an initially large 77™ can 
be compensated by an asymmetry in the other flavours with opposite sign. 
The most restrictive BBN bound on rj^^ applies then to the total asymmetry, 
an effect that can be seen graphically in figure [3] as a rotation of the allowed 
region from a quasi-horizontal one for zero mixing to an almost vertical region 
for sin^6'i3 = 0.04, in particular for the IH. In all cases depicted in figure [s] 
the allowed values of the asymmetries are mainly fixed by the ^He bound, 
which imposes that the value of 7),^^ at BBN must be very close to zero, while 
the data on primordial deuterium is crucial for closing the region. 

For values of 6*13 close to the upper limits set by experimental data, the 
combined effect of oscillations and collisions leads to an efficient mixing of 
all neutrino flavours before BBN. Therefore, the individual neutrino asym- 
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Figure 5: Largest values of iVcff from primordial neutrino asymmetries compatible with 
BBN at 95% C.L., as a function of 9i3 and the neutrino mass hierarchy. 

metries have similar values, approximately rj^^ ~ ri,y/3, and the BBN bound 
on the electron neutrino asymmetry applies to all flavours, and in turn to rj^y 
as considered in previous analyses [301 EB ESI ESj [SU |35l |36] . We find that 
for sin^^ig = 0.04 the allowed region at 95% C.L. is -0.17(-0.1) < r]^ < 
0.1(0.05) for neutrino masses following a normal (inverted) mass hierarchy. 
Note, however, that in the IH this result approximately holds for any value 
of sin^ within the favoured region by oscilation data, due to the resonant 
character of the conversions. Instead, as discussed in [23], in the NH even 
values of order \ri^\ ~ 0.6 are still compatible with BBN if ^13 is very small. 
The allowed regions of the total neutrino asymmetry are depicted in figure |4] 
as a function of the mixing angle ^13 and the mass hierarchy. 

Finally, the dependence of the largest value of A'^efr from neutrino asymme- 
tries in the region compatible with BBN, as a function of the neutrino mass 
hierarchy and the mixing angle 613 is reported in figure [5j If the true value 
of ^13 lies in the upper part of the region favoured by oscillation experiments 
(in particular T2K) or Amg^^ < 0, the presence of primordial asymmetries 
can not lead to a contribution to the radiation energy density Nes > 3.1 . 
On the other hand, for the NH and very small values of ^13, larger values of 
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Nes are still compatible with BBN data, up to 3.43 at 95% CL. 
4. Conclusions 

We have found the BBN constraints on the cosmological lepton number 
and its associated contribution to the radiation energy density, taking into 
account the effect of flavour neutrino oscillations. Once the other neutrino 
mixing parameters have been fixed by oscillation data, we have shown that 
pinpointing the value of ^13 is crucial to establish the degree of conversion of 
flavour neutrino asymmetries in the early Universe. 

We conclude that the most stringent BBN bound on the total neutrino 
asymmetry, \r]i,\ < 0.1, requires that reactor neutrino experiments in the 
near future [57], such as Double Chooz, Daya Bay or Reno, confirm that the 
value of the third neutrino mixing angle is such that sin^ ^13 > 0.03. This 
conclusion also applies to the whole allowed range of 613 values at 3(j [21 13] if 
neutrino masses follow an inverted hierarchy scheme. For smaller values of 
this mixing angle in the NH the BBN bound is relaxed, up to a maximum 
allowed range of —0.7 ^ rju ^ 0.6 for Oi^ = 0. 

Similarly, a measured value sin^ ^13 > 0.03 will imply that the maximum 
contribution of neutrino asymmetries to the radiation content of the Universe 
can not exceed A'^eff — 3.1, well below the expected sensitivity of the Planck 
satellite (0.4 at 2a) [3HI EHl SO], whose first data release on the anisotropies 
of the cosmic microwave background (CMB) is expected in one year or so. 

Finally, we emphasize that BBN remains the best way to constrain a 
potential cosmological lepton number, despite the present precision of the 
measurements of the CMB anisotropies and other late cosmological observ- 
ables. Bounds on the neutrino asymmetries with these data [H] H3| HI] 
do not improve those found in our work, but are of course sensitive to other 
neutrino properties such as their masses. 
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